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Learning Goals 


• Discuss some broad guidelines for turbulence models. 

• List of available turbulence models (calorically perfect gas) 

• Discuss the typical namelist parameters used. 

• Show some sections of fun3d.nml namelists used for turbulent flow 
simulations. 

• The detailed theory of turbulence models will not be covered in this 
session. 

• Pros and cons of each model will not be discussed either due to 
time limitations. 

- All of the models will likely work some of the time. 

- But none of the models will work all of the time. 
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The List 

Steady flow simulations 


• One-equation 

- Spalart-Allmaras (sa), Recherche Aerospatiale, No. 1, 1994. 

- Negative Spalart-Allmaras (sa-neg), ICCFD7-1902, 2012. 

• Two-equation 

- Menter-SST (sst),AIAAJ (32), 1994. 

- Menter-SST with vorticity source term (sst-v), NASA-TM-1 03975, 1992. 

- Menter-SST from 2003 (sst-2 0 03), Turbulence, Heat and Mass Transfer 4. 

- Wilcox k-omega (wiicox2006), AIAAJ (46), 2008. 

- Wilcox k-omega (wilcoxl 9 98), Turbulence Modeling for CFD, 1998. 

- Wilcox k-omega (wilcoxl988), AIAAJ (26), 1988. 

- Nonlinear k-omega (EASMko2 0 03-s),J Aircraft (38), 2001. 
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The List 

Steady flow simulations 


• Four-equation 

- Langtry-Menter transition model (gamma-ret-sst), AIAA-2005-0522. 

• Seven-equation 

- Wilcox Stress-omega RSM (wilcoxRSM-w2 0 0 6), Turbulence Modeling for CFD, 
2006. 

- SSGLRR-RSM (ssglrr-rsm-w2 012), AIM Journal, Vol. 53, No. 3, 2015, pp. 
739-755. 


Other references and detailed explanations of the models can be found 

at the turbulence modeling website: 



http : //turbmodels . larc . nasa . gov 
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The List 

Time accurate flow simulations 


• One-equation 

- Detached eddy simulations, (des, des-neg), TCFD (20), 2006. 

• Two-equation 

- Hybrid RANS-LES (hrles), AIAA-2008-3854. 
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http : 


General usage guidelines 

Do we even need to perform a turbulent flow simulation? 
Flow physics 

• What physics need to be simulated/predicted? 

• high speed flow -> possibly largely laminar 

• corner flow -> possibly anisotropic turbulence 

• blunt body wake -> possibly large eddy simulations 

Computational requirements 

• to evaluate the grid's resolution required for a certain accuracy 


physics 


accuracy 


resources 



Vwww. stanford.edu/ class/me 4 6 9b /handouts / turbulence . pdf , slide 51 


(fD i 
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CFD Process 


General usage guidelines 

• Appropriate spacing of the mesh on viscous solid walls must be 
used. 

- Generally accepted spacing is between .1 and 2.5 wall units. 

- Using wall functions, generally accepted spacing is between 0.1 and 250 wall 
units. 

- Many problems may have multiple scales, so no one physical distance for the 
first node spacing will suit the whole problem. 

• Generate a mesh with appropriate resolution to model the problem ( 
within the limits of the available computational resources ). 

- Try not to expand the mesh spacing too quickly away from a viscous wall. 

- Typically the more curvature in the physical geometry, the higher concentration of 
mesh. 

• One-equation models like Spalart-Allmaras tend to be very robust, 
cover a very wide range of flow situations and are a compromise 
between simplicity and accuracy. 

• Multi-equation models like the Menter-SST or RSM require more 
computational resources, but are more physically complete and can, 
possibly, add more accuracy to the solution... though YM MV. 
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General usage guidelines 


• Solutions to a steady state are adequate for many problems. 

• Depending upon the physics of the simulation, though, time- 
accurate solutions may be required. 
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Namelists 


fun3d . nml 


For turbulent flow simulations, depending upon the turbulence model 
and problem the following namelists within fun3d.nml are used. 

• &governing_equations 

• &turbulent_dif fusion_models 

• &spalart 

• &gammaretsst 
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Spalart-Allmaras 

fun3d . nml 


&governing_equations 

eqn_type = ’ cal_per_compress ’ 

viscous_terms = ’turbulent’ 

/ 

& turbul en t_di f f u s i on_mode 1 s 

turbulence_model = ' sa ' ! default 

! current 1-eqn options: ’sa-neg’, ’ des ’ , ’ des-neg’ 

turb_compr e s s_mode 1 = ’ none ’ 

! current options: ’ssz’ ! (Ref. AIAA-95-0863, Shur et al.) 

/ 
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Spalart-Allmaras 

fun3d . nml 


&spalart 
turbinf = 3.0 

! free stream value for spalart model 
ddes = . false . 

! for activating delayed DES model 
ddes_modl = .false. 

! Mod to DDES , Ref. AIAA Paper 2010-4001 
sarc = .false. 

! Ref. AIAAJ, Vol.38, No. 5, 2000, pp. 784-792. 
sarc_cr3 = 1.0 

! constant associated with SARC model 

/ 




N3D 

Fufy Uhstructured Navier-Stokes 


http://fun3d.larc.nasa.gov 


FUN3D Training Workshop 
June 20-21, 2015 


11 


Menter-SST 

fun3d . nml 


&governing_equations 

eqn_type = ’ cal_per_compress ’ 

viscous_terms = ’turbulent’ 

/ 

& turbul en t_di f f u s i on_mode 1 s 
turbulence_model = ’ sst ’ 

! other options: ’sst-v', ’sst-2003’, ’ gamma-ret-sst ’ 

! ’hrles’ 

/ 

&gammar e t s s t 

set_k_inf_w_turb_intsty_j?ercnt =0.2 ! (percent) 
set_w_inf_w_eddy viscosity =1.0 ! (nondim) 

transition_4eqn_on = . true . 

! toggles transition 

/ 
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Sample fun3d.nml 

Subsonic bump using S-A 

http : //turbmodels . larc . nasa . gov/bump . html 


fiproject 

pro ject_rootname = 'bump_31evelsdown_177x81 


/ 


&reference_physical_p roper ties 


ma ch_numbe r 
reynolds_number 
temperature 
tempera ture_units 

/ 


= 0.2 

= 3000000.0 
= 540.0 
= 'Rankine' 


&turbulent_dif fusion_models 
turbulence_model = ' sa ' 

/ 

&nonlinear_solver_parameters 


schedule_ 

_iteration 

= 1 

250 

schedule_ 

_cfl 

= 10. 

250. 

schedule 

cf lturb 

= 10. 

250. 


&boundary_condi tions 
total_pressure_ratio (3) = 1.02828 

total_temperature_ratio (3) = 1.008 
static_pressure_ratio (4) =1.0 

/ 
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Bump-in-channel Boundary Conditions, 



M=0.2, Re L = 3 million (L=1 ), T re) = 540 R 

4 

_ 

t 

symmetry 



- 

^ Pt/P rei =1 .02828, 



- 

Tt/T wl = 1.008, 




1 quantity from interior p/p _•( q, ► 


>- 
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2 
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Sample fun3d.nml 

Time accurate simulation using a S-A based DES model 


& turbulen t_di f f u s i onjnode Is 

turbulence model = ' des ' 


/ 


&nonlinear_solver_parameters 
time_accuracy = 

time_step_nondim = 

pseudo_time_stepping = 

subiterations = 

schedule_iteration = 

schedule_cfl = 

schedule cflturb = 


/ 


' 2ndorderOPT ' 
0.10 
' on' 

10 

1 100 
5 . 5 . 

5 . 5 


Details of running a time accurate simulations are covered in Session 1 1 . 
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EOF 

Turbulent flow simulations with Fun3D 


Several turbulence model options are available in VI 2.7 
Namelist nomenclature has been discussed. 

Caveats: 

Meshing and turbulence model decisions are highly 
dependent on the degree of fidelity and accuracy desired. 

The desired aspects, though, may not fit inside the resources 
available. 
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